Delocalization in harmonic chains with long-range correlated random masses 
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We study the nature of collective excitations in harmonic chains with masses exhibiting long- 
range correlated disorder with power spectrum proportional to l/fc", where k is the wave- vector 
of the modulations on the random masses landscape. Using a transfer matrix method and exact 
diagonalization, we compute the localization length and participation ratio of eigenmodes within the 
band of allowed energies. We find extended vibrational modes in the low-energy region for a > 1. 
In order to study the time evolution of an initially localized energy input, we calculate the second 
moment M2{t) of the energy spatial distribution. We show that M2{t), besides being dependent of 
the specific initial excitation and exhibiting an anomalous diffusion for weakly correlated disorder, 
assumes a ballistic spread in the regime a > 1 due to the presence of extended vibrational modes. 
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I. INTRODUCTION 

The role played by disorder on the nature of collec- 
tive excitations in condensed matter physics has been 
the subject of intensive studies due to its relevance 
in defining general transport characteristics [Q. Usu- 
ally, disorder induces localization of collective excita- 
tions thus degrading transport properties, an effect that 
is largely pronounced in low dimensions, fn particu- 
lar, the one-electron eigen-states in the one-dimensional 
Anderson model with site-diagonal uncorrelated disorder 
are exponentially localized for any degree of disorder . 
However, several one-dimensional models with correlated 
disorder have been proposed which exhibit delocalized 
states H, Recently, it has been shown that the 

one-dimensional Anderson model with long-range corre- 
lated disorder presents a phase of extended electronic 
states [H, 1^, ^ . These results have been confirmed by mi- 
crowave transmission spectra of single-mode waveguides 
with inserted correlated scatters ||] . 

The above results have motivated the study of further 
model systems that can be mapped onto the Anderson 
model and, therefore, expected to present a similar tran- 
sition between localized and extended collective excita- 
tions. Recently, a study concerning the one-dimensional 
quantum Heisenberg ferromagnet with exchange cou- 
plings exhibiting long-range correlated disorder reported 
some finite-size scaling evidences of the emergence of a 
phase of extended low-energy excitations . By using 
a renormalization group calculation the existence of such 
phase of extended spin- waves was confirmed and the scal- 
ing of the mobility edge with the degree of correlation was 
obtained It was also shown that, associated with the 
emergence of extended spin-waves in the low-energy re- 
gion, the wave-packet mean-square displacement exhibits 
a long-time ballistic behavior. 

The collective vibrational motion of one-dimensional 
disordered harmonic chains of N random masses can also 



be mapped onto an one-electron tight-binding model [|2j . 
In such a case, most of the normal vibrational modes 
are localized. However, there are a few low-frequency 
modes not localized, whose number is of the order of 
^/N , in which case the disordered chains behaves like the 
disorder-free system Futher, it was shown that 

correlations in the mass distribution produce a new set of 
non-scattered modes in this system Non-scattered 
modes have also been found in disordered harmonic chain 
with dimeric correlations in the spring constants [^ . 
By using analytical arguments, it was also demonstrated 
that the transport of energy in mass-disordered (uncorre- 
lated and correlated) harmonic chains is strongly depen- 
dent on non-scattered vibrational modes as well as on 
the initial excitation [ p6| . For impulse initial excitations, 
uncorrelated random chains have a superdiffusive behav- 
ior for the second moment of the energy distribution 
[M2{t) oc t^'^], while for initial displacement excitations 
a subdiffusive spread takes place [M2{t) cx t^-^\. The 
dependence of the second moment s pre ad on the initial 
excitation was also obtained in Ref. |]1^. Moreover, cor- 
relations induced by thermal annealing have been shown 
to enhance the localization length of vibrational modes, 
although they still present an exponential decay for dis- 
tances larger than the thermal correlation length [ p^ . 
Recently the thermal conductivity on harmonic and an- 
harmonic chains of uncorrelated random masses |p^ , as 
well as of the chain of hard-point particles of alternate 
masses [ pO| , has been numerically investigated in detail. 
The main issue here is whether the systems display fi- 
nite thermal conductivity in the thermodynamic limit, a 
question that remains controversial |2|] . 

In this paper we extend the study of collective modes in 
the presence of long-range correlated disorder for the case 
of vibrational modes. We will consider harmonic chains 
with long-range correlated random masses assumed to 
have spectral power density S oc By using a trans- 

fer matrix calculation, we obtain accurate estimates for 



2 



the Lyapunov exponent, defined as the inverse of the de- 
gree of locahzation of the vibrational modes. We show 
that, for a > 1, this model also presents a phase of ex- 
tended modes in the low frequency region. This result 
is confirmed by participation ratio measurements from 
an exact diagonalization procedure and finite size scaling 
arguments. The spatial evolution of an initially localized 
excitation is also studied by computing the spread of the 
second moment of the energy distribution, M2(t). We 
find that, associated with the emergence of a phase of 
delocalized modes, a ballistic energy spread takes place. 



II. FORMALISM 

We consider a disordered harmonic chain of N masses, 
for which the equation of motion for the displacement u„ 
of the n-th mass with vibrational frequency u is JTb] , Q 



(1) 



Here m„ is the mass at the n-th site and /?„ is the spring 
constant that couples the masses to„ and run+i- We 
use units in which /?„ = !. In the present harmonic 
chain model, we take the masses run following a random 
sequence describing the trace of a fractional Brownian 
motion g || ||: 
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where k is the wave- vector of the modulations on the ran- 
dom mass landscape and cjjk are N/2 random phases uni- 
formly distributed in the interval [0, 27r]. The exponent a 
is directly related to the Hurst exponent H {a = 2H + 1) 
of the rescaled range analysis. In order to avoid vanish- 
ing masses we shift and normalize all masses generated 
by Eq. (2) such to have average value < m„ >= 5 and 
variance independent of the chain size (< Am„ >= 1). 
Using the matrix formalism, Eq. (1) can be rewritten 

as 
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For a specific frequency w, a 2 x 2 transfer matrix T„ 
connects the displacements at the sites n — 1 and n to 
those at the site n -I- 1: 
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FIG. 1; Normalized density of states (DOS) as a function 
of Lo'^ obtained using Dean's method. The chain length is 
— for all cases. The DOS becomes less rough as a is 
increased. For a = 1.5 it displays a non-fluctuating part near 
the bottom of the band. 
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FIG. 2: Lyapunov coefficient 7 as a function of lj for a = 
(uncorrelated random chain) and N = 2x 10^ sites. The Lya- 
punov coefficient is finite for non-zero frequencies (localized 
states) and vanishes as 7 oc cj^, — * (inset). 



The localization length of each vibrational mode is taken 
as the inverse of the Lyapunov exponent 7 defined by p^ , 

MM 
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Once the initial values for uq and ui are known, the value 
of Un can be obtained by repeated iterations along the 
chain, as described by the product of transfer matrices 
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where c(0) — ("^) is a generic initial condition. Typically, 
2 X 10^ matrix products were used to calculate the Lya- 
punov exponents. The nature of the vibrational modes 
can also be investigated by computing the participation 
ratio ^, since it displays a dependence on the chain size 
for extended states and is finite for exponentially local- 
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ized ones. ^ is defined by ^ 

= ^r^, (7) 

where the displacements it„ are those associated with 
an eigenmodes a; of a chain of N masses and are ob- 
tained by direct diagonahzation of the N x N secular 
matrix A defined by A^^i = + (3i_i)/mi = 2/mi, 
A^^i+i = Ai+i^i = /^(m^mj+i) = l/{mimi+i), and all 
other Aij = The participation ratio calcula- 

tions were averaged over 100 samples. 

We compute the Lyapunov exponents and the partic- 
ipation ratio for several values of the correlation expo- 
nent a, and obtain the density of states (DOS) using the 
numerical Dean's method ||2^. Strong fluctuations in 
the DOS are related to the presence of localized states, 
whereas smooth a DOS is usually connected with the 
emergence of delocalized states Q . In Fig. 1 we show 
the normalized DOS for chains with N — 10^ sites, and 
notice that it becomes less rough as a is increased. In 
Fig. 2 we display the plot of 7 versus for a — (un- 
correlated random masses). The Lyapunov coefficient 
is finite for all frequencies and vanishes at = as 
7 (X w^, in agreement with Ref. ||T^ (see inset of Fig. 2). 
The scaled participation ratio ^/N as a function of is 
shown in Fig. 3: ^(cj^ = 0)/N remains finite in the ther- 
modynamic limit, whereas for any finite frequency the vi- 
brational modes are localized with ^/TV^OasiV^cx), 
in agreement with the above results obtained from the 
Lyapunov coefficient calculations. To investigate the ef- 
fect of weak long-range correlated disorder, we present 
in Fig. 4(a) the Lyapunov coefficient as a function of lj^ 
for a — 0.75 and = 2 x 10^. In spite of 7 being very 
small in the bottom of the band, the scaled participation 
ratio for lu > vanishes in the thermodynamic limit [see 
Fig. 4(b)]. Therefore, all modes with a; > are still lo- 
calized, a feature that holds for any < a < 1. However, 
the nature of the low-frequency modes changes qualita- 
tively for a > 1. In Fig. 5(a) we show 7 versus for 
a = 1.5 and N = 2x 10^ sites. The Lyapunov coefficient 
vanishes within a finite range of frequency values, thus re- 
vealing the presence of extended vibrational modes. The 
scaled participation ratio ^/N [see Fig. 5(b)] displays a 
well defined data collapse, confirming that the phase of 
extended low-frequency vibrational modes is stable in the 
thermodynamic limit. 



III. ENERGY TRANSPORT 

In order to study the time evolution of a localized en- 
ergy pulse, we calculate the second moment of the energy 
distribution [^ . This quantity is related to the ther- 
mal conductivity by Kubo's formula [|l^, The classi- 
cal Hamiltonian H for an harmonic chain can be written 
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FIG. 3: Scaled participation ratio £_/N as a function of oj^ 
for a — (uncorrelated disorder). From top to bottom. A'' = 
200, 400, 800, 1600. For vibrational modes with w > 0, ^/N 
as N diverges, thus confirming their localized nature. 
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FIG. 4: (a) Lyapunov coefficient 7 versus lj^ for a = 0.75 
and N = 2 X lO'' sites, (b) Scaled participation ratio £,/N 
as a function of oj^ for a = 0.75. From top to bottom, 
iV = 200, 400, 800, 1600. In spite of 7 being very small in 
the bottom of the band, all modes with a; > are localized. 

as 

H = J2hnit), (8) 
where the energy /i„(i) at the site n is given by 

hnit) = ^ + ^[iQn+l-Qnf + iQn-Qn-lf] ■ (9) 

Here P„ and Qn define the momentum and displacement 
of the mass at the n-th site. The Hamilton's equations 
are 

Pn(t) = -— - = /3„[(Q„+i - Qn) - [Qn ~ Qn-l)] (10) 
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FIG. 5: (a) Lyapunov coefficient 7 versus lu'^ for a — 1.5 and 
N = 2 X 10^ sites. Tlie Lyapunov coefficient vanisfies witfiin 
a finite range of frequency values, tfius revealing the presence 
of extended vibrational modes, (b) Scaled participation ratio 
S^/N as a function of uj'^ for a = 1.5. From top to bottom, 
TV = 200, 400, 800, 1600. The phase of extended vibrational 
modes is confirmed by the size independent plateau in the 
low- frequency region. 



and 
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The fraction of the total energy H at the site n is 
given by hn{t)/H and the second moment of the energy 
distribution, M2{t), is defined by p6| 



N 



(12) 



n=l 



where an initial excitation is introduced at the site tt-q at 
t = Q . Using the fourth-order Runge-Kutta method, we 
solve the differential equations for P„(i) and Qn{t) and 
calculate M2{t). The second moment of the energy dis- 
tribution M2{t) has the same status of the mean-square 
displacement of the wavepacket of an electron in a crys- 
tal p^ . In harmonic chains with an initial impulse ex- 
citation, the energy spread is faster than that in chains 
with an initial displacement excitation |l^, We cal- 
culate M2 (t) for several a values and two kinds of initial 
excitation: impulse excitation and displacement excita- 
tion. 
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FIG. 6: Scaled energy second moment M2{t)/t^'^ versus time 
t for a = (dotted line) and a — 0.75 (dashed line) with 
initial impulse excitation. For < a < 1 only superdiffusive 
behavior is found for long times. 



0.75 (dashed line). These results have been obtained af- 
ter an initial impulse excitation, Pnoi^ = 0) = 5n„ N/2- 
In our calculations for a = 0, the self-expanded chain 
method with initial chain size N = 1000 was used to 
minimize end effects. Throghout the numerical integra- 
tion process we kept the fraction of the total energy H 
at the ends of the chain [ho{t)/H and hN{t)/H] smaller 
than 10"^°'' for all times. As shown in Fig. 6, we find a 
long-time superdiffusive behavior for a = 0, in agreement 
with previous analytical and numerical results for energy 
transport in harmonic chains with uncorrelated random 
masses under an impulse initial excitation |l^ . In con- 
trast, for a > we cannot use the self-expanded chain 
method due to the long-range character of the mass cor- 
relations. Therefore, chains with N = 10000 masses were 
considered, and the runs stopped whenever the fraction 
of the total energy at the chain ends achieved 10"^°". 
For a = 0.75 the time-dependence of the scaled energy 
second moment M2{t)/t^'^ typically represents a weak 
long-range correlated case. In such a case we also find su- 
perdiffusive behavior for long times. On the other hand, 
in the strong correlated regime, a > 1, a breakdown in 
the superdiffusive behavior sets up. Fig. 7 shows the 
time-dependence of the scaled energy second moment, 
M2(^)/^^ for a = 1.5 (dotted line) and a = 2.0 (dashed 
line). Associated with the emergence of extended vibra- 
tional modes in the low-energy region, the second mo- 
ment M2{t) displays a long-time ballistic behavior. 



B. Displacement excitation 



A. Impulse excitation 

In Fig. 6 we present the scaled second moment 
M2{t)/t^'^ versus time <: for a = (dotted line), which 
corresponds to the uncorrelated random chain, and a — 



The long time behavior of the second moment M2 (t) in 
uncorrelated random chains with initial displacement ex- 
citation is significantly different from the corresponding 
behavior with impulse initial excitation p^ , . Analyt- 
ical calculations predict that M2(t) oc a result that 
has been corroborated by numerical techniques p6| . For 
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FIG. 7: Scaled second moment M2{t)/t^ versus time t for a = 
1.5 (dotted line) and a = 2.0 (dashed line) with initial impulse 
excitation. Results were obtained by numerical integration in 
chains with A'^ = 10000 sites. A ballistic behavior is found 
after an initial transient. 
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FIG. 8: The scaled second moment M2{t)/t^'^ versus time 
t for a = O(dotted line) and a = 0.75 (dashed line) with 
initially given displacement excitation. The subdiffusive be- 
havior is found for long times. 



a = Q wc indeed reproduce this behavior, as shown in 
Fig. 8 for the scaled second moment M2{t)/t^'^ versus 
time t with initial displacement excitation Qn^ (t = 0) = 
^no,N/2- Again, we find that this asymptotic subdiffu- 
sive behavior remains true for < a < 1 (dashed hne in 
Fig. 8). For strong correlations {a > 1), which induce the 
emergence of new extended vibrational modes in the low- 
energy region, the energy transport is faster than in the 
subdiffusive regime and again assumes a ballistic nature, 
as shown in Fig. 9. 



IV. SUMMARY AND CONCLUSIONS 
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FIG. 9: The scaled second moment M2{t)/t^ versus time t for 
a = 1.5 (dotted line) and a — 2.0 (dashed line) with initially 
given displacement excitation. Results obtained by numerical 
integration of diferential equation for chain with A'^ = 10000 
sites. A ballistic behavior is found for all times. 



and exact diagonalization, we have computed the local- 
ization length and the participation ratio of all normal 
modes. Our results indicate that in the strong correla- 
tions regime, a > 1, there is a phase of extended low- 
energy vibrational modes. In this sense, long-range cor- 
relations in the mass distribution induce the emergence of 
a delocalization transition in harmonic chains similar to 
the one observed to occur with one-magnon excitations 
in ferromagnetic chains with random couplings ^] 
and with one-electron eigen-states in the random hopping 
Anderson model j7j. We have also studied the energy 
transport in this harmonic chain model. The spread of 
the energy second moment M2 (t) is shown to be strongly 
dependent on the existence of non-scattered vibrational 
modes and initial excitation. We have also found that, as- 
sociated with the emergence of a phase of low-energy ex- 
tended collective excitations, M2(t) displays a crossover 
from an anomalous sub- or super-diffusive regime (de- 
pending on the initial impulse or displacement excita- 
tion, respectively) to an asymptotic ballistic behavior. 
The above findings indicate that the thermal conductiv- 
ity can be strongly influenced by the presence of long- 
range correlations in the random distribution of masses 
and we hope that the present work will stimulate further 
studies along this direction. 
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In this paper we have studied the nature of collective 
vibrational modes in harmonic chains with long-range 
correlated random masses m„, with spectral power den- 
sity S oc l/k". By using a transfer matrix method 
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